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1.  INTRODUCTION 

One  who  seeks  to  familiarise  himself  with 
adaptive  antenna  systems  for  interference  rejection 
is  confronted  by  a  vast  literature  containing  many 
different  proposed  systems.  However,  there  are 
only  a  few  fundamentally  different  concepts  under¬ 
lying  the  variations.  In  this  paper,  these  concepts 
and  their  relationships  are  stressed.  Adaptive  sys¬ 
tems  that  are  not  intended  for  interference  rejec¬ 
tion,  such  as  retrodirective  arrays,  are  not  con¬ 
sidered. 

Special  attention  is  given  to  the  sidelobe  can¬ 
celler,  the  archetypical  adaptive  antenna  system. 
The  initial  mathematical  analysis  is  heuristic  in 
order  to  clarify  the  basic  principles  of  operation.  It 
is  shown  that  the  adaptation  in  a  sidelobe  canceller 
can  be  interpreted  not  only  as  noise  cancellation, 
but  also  as  adaptive  beam  forming  and  null  steer¬ 
ing.  The  analysis  of  the  sidelobe  canceller  leads 
naturally  to  consideration  of  the  multiple  sidelobe 
canceller.  A  slight  modification  of  the  sidelobe 
canceller  yields  a  notch  filter. 

After  the  initial  heuristic  analysis  of  the  side¬ 
lobe  canceller,  the  remainder  of  this  paper  contains 
a  more  rigorous  mathematical  analysis.  Alterna¬ 
tive  adaptive  systems  that  can  handle  wide-band- 
width  signals  are  developed  from  the  classical 
theory.  There  are  two  forms  of  the  adaptive  ele¬ 
ments.  In  one  form,  the  weight-adjustment  mechan¬ 
ism  responds  to  the  system  output;  in  the  other 
form,  the  mechanism  responds  to  the  difference 
between  the  output  and  the  desired  response.  The 
inputs  to  the  adaptive  elements  may  be  derived 
from  tapped-delay-line  or  frequency-domain  array- 
procesiiing  filters. 

The  detailed  design  of  an  adaptive  antenna 
system  is  dependent  upon  the  performance  cri¬ 
terion  selected.  Two  of  the  most  useful  criteria  are 
the  signal-to-noise  ratio  and  the  mean-square- 
error  criteria.  The  latter  criterion  combined  with 
the  method  of  steepest  descent  leads  to  the  Widrow- 
Hoff  algorithm.  This  algorithm  is  often  imple¬ 
mented  in  practical  systems. 

The  adaptive  noise  canceller  differs  from  classi¬ 
cal  adaptive  systems  primarily  in  that  the  estimate 


of  the  desired  response  is  derived  from  a  separate 
antenna,  called  the  primary  antenna,  rather  than 
from  an  internt J  generator.  The  sidelobe  canceller 
is  a  special  tyfie  of  adaptive  noise  canceller.  The 
adaptive  notch  filter,  which  removes  periodic 
interference  s'gnals  from  wide-bandwidth  desired 
signals,  is  another  variation  of  the  adaptive  noise 
canceller. 

Classical  adaptive  systems  and  the  adaptive 
noise  canceller  sometimes  exhibit  unintentional 
cancellation  of  the  desired  signal.  To  preclude  such 
a  possibility,  the  constrained  minimum  power 
criterion  can  be  used  in  the  design  of  an  adaptive 
system.  The  resulting  system  automatically  limits 
the  cancellation  of  the  desired  signal,  while  still 
adaptively  filtering  the  interference. 

2.  SIDELOBE  CANCELLER 

An  adaptive  antenna  system  automatically 
monitors  its  output  and  adjusts  its  parameters 
accordingly.  It  does  so  in  order  to  reduce  the 
impact  of  interference  that  enters  through  the 
sidelobes,  or  possibly  the  mainlobe,  of  its  antenna 
radiation  pattern,  while  still  allowing  reception  of 
an  intended  transmission.  The  design  of  an  adap¬ 
tive  antenna  system  requires  little  a  priori  knowl¬ 
edge  of  the  signal  or  interference  characteristics. 

The  sidelobe  canceller  is  a  classic  example  of  an 
adaptive  antenna  system.  It  not  only  is  of  practical 
importance,  but  also  provides  an  introduction  to 
the  fundamental  concepts  of  adaptive  antenna 
systems. 

Figure  1  shows  a  version  of  a  sidelobe  canceller. 
The  primary  and  reference  signals  are  outputs  of 
two  separate  antennas  or  two  different  groups  of 
elements  in  a  single  phased-array  antenna.  It  is 
intended  that  the  reference  signal  provide  an 
estimate  of  the  interference  in  the  primary  signal. 
After  suitable  processing,  this  estimate  is  sub¬ 
tracted  from  the  primary  signal.  As  a  result,  the 
interference,  which  may  have  entered  through  the 
sidelobes  of  the  primary  antenna,  is  reduced  or 
eliminated  by  cancellation. 

Ideally,  the  reference  signal,  X,(t),  has  a  large 
interference  component  and  a  small  desired  signal 
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component,  whereas  \0(/)  may  have  a  much  larger 
icsired  signal  component.  The  quart er*wavelength 
delay  shown  in  figure  1  produces  a  signal,  Xj(t), 
that  is  in  phase  quadrature  with  X,(f)  (a  quad¬ 
rature  hybrid  could  be  used  instead  of  the  delay). 
The  weight  functions.  W,(t)  and  W,(<)«  regulate  the 
amounts  of  X,(<)  and  X,(t)  that  are  subtracted  from 
X„(f).  The  relative  magnitudes  of  the  weight  func¬ 
tions  determine  the  magnitude  and  phase  of  the 
total  waveform  that  is  subtracted  from  X#(/l  If  the 
magnitudes  and  phases  of  the  interferer  com¬ 
ponents  of  X,(t)  and  the  total  subtracted  waveform 
are  equal,  the  interference  is  cancelled  and  does  not 
appear  in  the  output.  y(t).  If  the  cancellation  is 
nearly  comp’e.e,  the  weight  functions  are  nearly 
constants;  if  it  is  not.  the  weight  functions  vary  in 
such  a  w^y  that  the  total  subtracted  waveform 
gradually  becomes  a  facsimile  of  the  interference 
component  of  X„(l). 

In  the  following  analysis,  we  make  no  assump¬ 
tions  about  the  directionality  of  the  primary  an¬ 
tenna  beam  or  the  reference  antenna  beam.  How¬ 
ever,  it  is  highly  desirable  to  have  the  primary 
beam  point  in  the  direction  of  the  desired  signal. 
The  beam-steering  network  forms  a  beam  in  the 
appropriate  direction  by  using  various  types  of  a 
priori  information.  If  the  antenna  is  part  of  a  radar 
system,  the  information  may  be  radar  return 
characteristics.  If  the  antenna  is  part  of  a  commun¬ 
ication  system,  the  information  may  be  the  char¬ 
acteristics  of  a  pilot  signal  transmitted  along  with 
the  message  signal.  The  antenna  system  then  locks 
onto  the  pilot  to  form  a  beam  in  the  direction  of  the 
transmitter. 


2.1  Steady-State  Operation 

In  general,  the  desired  signal  received  by  the 
primary  antenna  has  the  form 

s(f)  =  As(f)  cos  [w„f  +  0,(0]  ,  (1) 

where  aj0  is  the  carrier  frequency ,  and  A,(f)  and 
0,(0  are  modulation  functions.  We  assume  that 
interference  (which  may  he  due  to  jamming),  after 
passage  through  the  bandpass  filter  of  the  primary 
branch,  has  the  form 


J(0  =  Ay(0  cos  [co,l  +  0j(O]  , 

when*  to,  is  the  carrier  frequency,  and  Ay(l)  and 
0y(O  an*  modulation  functions.  Neglecting  the 
thermal  noise,  the  output  of  the  bandpass  filter  of 
the  primary  branch  is 

X#(0  =  A,(0  cos  l«V  +  <A*(0]  + 

Ay(0  cos  [w,r  +  0y(O]  •  (3) 

Because  of  a  possible  difference  in  the  radiation 
patterns  between  the  two  antennas,  the  signals 
received  by  the  reference  antenna  may  experience 
different  amplifications  than  the  same  signals  do 
when  received  by  the  primary  antenna.  Each  signal 
has  a  different  arrival  time  and,  hence,  a  phase  shift 
at  the  reference  antenna  relative  to  the  same  signal 
at  the  primary  antenna.  We  assume  that  the 
antennas  are  close  enough  that  the  difference  in 
arrival  time  of  a  signal  at  the  two  antennas  is  much 
less  than  the  inverse  signal  bandwidth.  Thus,  the 
modulation  functions  are  negligibly  uffected  by 
this  difference.  Vi  e  conclude  that  the  output  of  the 
bandpass  filter  of  the  reference  branch  can  be 
represented  bv 

X,(t)  =  C,A»(r)  cos  [to„f  +  0s(r)  +  0,(r)]  + 

CjAy(t)  COS  [(0,f  +  0y(f )  +  0>(t)]  •  (D 

where  0,  and  ft.  are  phase  angles,  and  C,  and  (\  are 
real  constants.  If  the  sources  of  the  desired  signal 
and  the  interference  are  separated  geometrically, 
then  0,  5^  03.  If  the  sources  are  n  obile,  then  0,  and 
t?2  are  functions  of  time.  The  magnitudes  of  the 
primary  antenna  gains  in  the  directions  of  the 
desired  signal  and  the  interference  are  denoted  by 
CpS  and  l wpj.  respectively.  The  magnitudes  of  the 
reference  antenna  gains  in  the  directions  of  the 
desired  signal  and  the  interference  are  denoted  by 
C.rs  and  (>,y.  respectively.  From  these  definitions, 
we  have 


The  integrators  are  designed  to  integrate  over 
the  time  interval  I  =  |f  —  T.f],  where  T  is  such 
that 
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w,T  »  1  .  (6*) 

We  utumr  that  the  bandpass  filter  is  narrow¬ 
band  so  titat 

»  2rrB  .  (6b) 

where  B  is  the  bandwidth  of  the  bandpass  filters  in 
herti.  Thus,  the  bandwidth  due  to  the  message 
modulation  is  much  less  than  the  carrier  fre¬ 
quency.  For  practical  values  of  ay  the  bandwidths 
associated  with  the  time  variations  of  0,  and  0, 
are  also  much  less  than  uy 

The  quarter- wave  length  delay  is  designed  to 
introduce  a  90-degrev  phase  shift  in  the  intended 
transmission.  Since  w,  -  u,  <  2rrB.  equation 
(6b)  indicates  that  this  delay  also  introduces  nearly 
a  90-degree  phase  shift  in  the  interference,  al¬ 
though  it  does  not  significantly  affect  the  modu¬ 
lation  waveforms.  Thus, 

X,(r)  -  C,A,(/)  sin  [ay  4-  0,(r)  +  0,(t)]  + 

CjAy(f )  sin  [w,t  +  0y(/)  +  0j(/)]  .  (?) 

The  signal-to-interference  ratio  at  the  primary 
input  is  defined  by 


tJ 

f,2(u)du 

fj  ( u)du 

f  A,1  du  +  ft  A*  cos(2woU  +  2 0,)du 
— - - -  (») 

Jl  Ay  du  +  Ji  Ay  cos(2u>,<i  +  2  <t>j)du 

In  general,  p,  is  a  (unction  of  time.  The  steady  state 
is  defined  to  exist  when  this  ratio  is  nearly  constant 
over  time  intervals  of  duration  T  or  longer. 

W  hen  the  steady  state  is  reached,  equations  (6) 
imply  that  the  second  terms  in  the  numerator  and 
denominator  are  negligible  compared  to  the  first 
terms.  Thus, 


where  p*  is  the  ratio  that  would  exist  if  the  primary 
pattern  were  omnidirectional.  Similarly,  the  signal- 
to-interference  ratio  at  the  reference  input  is  found 
to  be 


Although  p,,  pn  and  p„  are  functions  of  time  in 
general,  they  are  nearly  constant- during  the  steady 
state. 

From  figure  1,  the  output,  y(t),  and  the  weight¬ 
ing  functions.  W',(t)  and  W  2(t),  are  given  by 

y(t)  -  X,(r)  -  W,(i)X,(i)  -  WrfOXM  . 

(ID 

w,(l)  -  tiff  y(u)X,(u)du  ,  (12) 

W,(t)  =  K  f  ylu)X,(u)du  ,  (13) 

where  the  constant  K.  is  th*  gain  of  a  linear 
amplifier.  Substituting  equation  (11)  into  equation 
(12),  we  obtain 

W,(f)  -K.fi  Xt(u)X,(u)du  - 

K  fi  W/,(u)X,2(u)du  ~ 

K./i  W|(u)X,(u)Xj(u)rfii  .  (14) 

We  shall  show  that  W ,(/)  and  W  2(t)  are  nearly 
constants  over  interval  I  if  pn  9 and  6}  are  nearly 
constants  over  interval  I.  In  this  case,W\(/)  »nd 
Wj(r)  can  be  removed  outside  the  integrals  with 
negligible  error,  so  that 

W,(r)  =  Kfx^X^du- 

KW,(i)  f  X’(u)du  ~ 

KW'j(f)  f  X,|u)Xj(u)du  .  (15) 
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In  a  annular  manner,  we  obuin 

W,(t)  -  K  /,  Xt(u)Xt(u)du  -  KW,(t) /,  X}( u)du 
-  KW,(i) /,  Xt(u)X,(u)du  .  (16) 

Equations  (6)  and  simple  trigonometry  show  that 


w,(l)  *  i 

/,  X«(u)X,(u)du 

(24) 

fl  X*  (u)du 

w,(t)  -  - 

fl  Xju)X,^)du 

(25) 

fl  X|  (u)du 

/,X5(u)du  *  /,  X?  (u)du  (17) 

during  the  steady  state.  Using  this  approximation 
to  solve  equations  (IS)  and  (16)  simultaneously 
yields 


Thus,  the  weighting  functions  are  essentially 
normalised  cross-correlation  functions.  From 
equations  (3),  (4),  (6),  and  (7),  we  obtain 

f\  X\(u)du  «  f  |  A,(u)  +  2,J^U)  + 

C|CjA,(«)Aj(u) 

COs[(w,  -  W,)tt  +  0,(li)  -  0y(u)  +  0,  -  >J  1  du  , 

'  (26) 


In  the  above  we  define 

v,  *  Y  +/tx»<'“  * 

(20) 

V,  =  /,  X,(u)Xj(u)du  , 

(21) 

V,  =  /,  X.(u)X,(u)du  . 

(22) 

V«  =  f  X,(u)X,(u)du  . 

(23) 

Evaluating  the  integrals  with  equations  (3),  (4), 
and  (7),  and  substituting  equations  (18)  and  (19) 
into  equation  (11)  gives  an  expression  for  the 
output  in  terms  of  the  desired  signal  and  inter¬ 
ference. 

In  most  cases  of  interest,  some  simplification  is 
possible.  If  T  and  K  are  sufficiently  large,  the  first 
term  on  the  right-hand  side  of  equation  (20)  can  be 
ignored.  Substituting  equations  (4)  and  (7)  into 
equation  (21)  shows  that  V,  is  very  small  during 
steady-state  operation.  Consequently,  it  is  rea¬ 
sonable  to  neglect  V,  in  equations  (18)  and  (19). 
As  a  result  of  these  approximations,  we  have 


fl  X,(u)X,(u)du  =  /|2‘  AJ(u)  cos  0,  + 

A/(u)  cos  0,  +  2‘  A,(u)A,<u) 
cos  [(w.  -  W,)u  +  0,(u)  -  0y<u)  +  0,]  + 

~  K(u)Ajiu) 

COSTCO,  -  W,)u  -)•  0,(u)  -  0,(u)  —  0,J  |  du  , 

(27) 

fl  X,(u)X,(u)rfu  =  Ai(“)  sin  + 

^  Ay(u)  sir  0,  +  A,(u)Ay(u) 
sin[(w0  -  w,)u  -  0,(u)  -  0y(u)  +  0|]  ~ 
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^2Af(u)A,<ii) 

8*n((w*  “  +  0/u)  “  0/«)  ~  0i]  |  du 

1  (28) 

In  general,  since  the  desired  signal  is  unsyn- 
chronixed  with  the  interference,  we  expect  that 
during  stef  -state  operation, 

J,(u.  0, )  J,(u.  0,)  du  =*  0  .  (29) 

where  a/u,  0,)  is  the  desired  signal  with  an  arbi¬ 
trary  phase  shift  of  0,.  and  J,(u.  0d  is  the  interfer¬ 
ence  signal  with  an  arbitrary  phase  shift  of  0,.  If 
equation  (29)  is  valid,  the  desired  signal  and 
interference  are  said  to  be  uncorrelated  in  the  time- 
average  sense. 


PrA/t)  cos  [«,*  +  0/1)]  - 
Or  a  A/l)  eoa  [ci»,t  +  0/1)  -  0,  +  0,]  (32) 

where 


The  signal-to-interference  ratio  at  the  output  of 
the  sidelobe  canceller,  p*  is  calculated  as  the  ratio 
of  the  average  power  over  interval  I  in  the  first  two 
terms  to  the  average  power  over  interval  I  in  the 
last  two  terms  of  the  right  side  of  equation  (32).  If 
y(t)  s*  0,  the  result  is  the  remarkably  simple 
formula, 

P*  *  —  •  (34) 

Pr 


Equation  (29)  implies  that  only  the  first  two 
terms  of  'quatiu.is  (26)  through  (28)  need  be 
retained.  Using  equations  (9)  and  (10),  equations 
(24)  through  (28)  yield 


W,  = 


C,  COS  6}  +  CjPrCOS  Bt 

C.C/l  +  pr ) 


(30) 


W,  = 


C|  sin  0,  +  Cj  pr  sin  0, 
C,Ci(l  4-  pr) 


(31) 


These  equations  indicate  that  W,  and  W}  are  ap¬ 
proximately  constant  over  interval  I  if  pn  0„  and  0, 
are  approximately  constant  over  interval  I .  Thus, 
the  initial  assumption  that  the  weighting  functions 
are  constant  has  been  verified  as  consistent  with 
the  other  approximations  made  in  deriving  these 
equations. 


We  now  substitute  into  equation  (11)  and 
employ  trigonometric  identities.  Using  equations 
(9)  and  (10),  the  final  result  is 

(1  +  Or)  y(t)  =  A/*)  cos [to/  +  0/1)]  - 
cF  A/I)  cos  [<v  +  0/1)  +  0,  -  0,]  + 


Thus,  the  interference  component  of  the  primary 
signal  has  been  nearly  cancelled  if  pr  <  <  1.  We 
conclude  that  the  output  signal  distortion  is  small 
when  the  signal  power  at  the  reference  antenna  is 
relatively  low.  An  important  approximation  made 
in  deriving  equations  (32)  and  (34)  is  that  the 
thermal  noise  is  negligible. 

If  the  interference  source  is  almost  directly 
behind  the  source  of  the  desired  signal,  then 
a  —  1  and  0,  —  02.  Since  equation  (32)  indi¬ 
cates  that  y(t)  -  0,  the  output  signal  is  buried  in 
the  thermal  noise.  We  conclude  that  the  sidelobe 
canceller  is  ineffective  in  this  case. 

A  strong  desired  signal  in  the  reference  channel 
can  induce  cancellation  of  the  desired  signal  at  the 
output.  Specifically,  if  p,  i  pr  >  1,  equation  (34) 
gives  p0  <  pj.  Consequently,  the  adaptive  system 
causes  a  performance  degradation  relative  to  the 
j»erformance  of  the  primary  antenna  operating 
alone.  Equation  (10)  indicates  that  designing  the 
system  so  that  Gra  «  is  helpful  in  ensuring 
small  values  of  pr.  Thus,  a  reference  antenna  with  a 
beam  that  points  in  the  direction  of  the  interference 
is  highly  desirable.  One  possible  implementation  is 
to  use  a  beam-steering  network  that  enables  a 
reference  beam  to  search  for  interfering  signals  of 
high  power  or  special  characteristics. 
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Sine*  the  dirretiun  of  the  interfe rence  is  usually 
unknown  a  priori,  and  since  there  may  be  two  or 
more  spatially  separated  interference  sources  pres¬ 
ent,  Gn  «  G,,  might  be  achieved  if  we  deploy 
several  reference  antennas,  each  with  a  directional 
beam,  as  shown  schematically  in  figure  2.  The 
reference  beams  point  in  various  directions  rela- 


small  that  pr  >  >  1 .  To  prevent  this  disaster,  the 
reference  input  can  be  monitored  and  its  propaga¬ 
tion  blocked  unless  its  amplitude  or  power  exceeds 
a  fixed  threshold .  Thus,  the  adaptive  mechanism  is 
disabled  unless  the  interference  is  strong  enough  to 
warrant  its  use. 


Figure  2.  Multiple  sidelobe  canceller. 


live  to  the  primary  beam.  The  output  of  each 
reference  antenna  feeds  a  cancellation  network 
similar  to  that  of  figure  1 .  A  configuration  of  this 
type  is  known  as  a  multiple  sidelobe  canceller  since 
it  provides  the  possibility  of  cancelling  multiple 
interfering  signals  entering  at  two  or  more  different 
points  of  the  sidelobes  of  the  primary  pattern.  A 
less  expensive  multiple  sidelobe  canceller  uses 
omnidirectional  reference  antennas.1 

Even  if  G,s  «  Grj,  serious  cancellation  of  the 
desired  signal  can  occur  if  the  interference  is  so 
- 

ft  /).  Strinbrrt,  Primaries  of  iperuinr  unW  4rr«v  Svslrm  m. 


If  the  reference  antenna  is  disconnected  and  an 
unmodulated  carrier  is  injected  as  X|(t),  the  adap¬ 
tive  antenna  system  behaves  as  a  notch  filter  at  the 
carrier  frequency.  To  demonstrate  this  effect,  we 
set  C,  =0  and  assume  that  Aj(t)  and  are 
constants  in  equations  (2)  and  (4).  If  equation  (29) 
is  valid,  we  obtain  y(f)  =  s(f),  indicating  com¬ 
plete  cancellation  of  the  interfering  carrier.  How¬ 
ever,  equation  (29)  is  not  valid  if  As(t)  and  0.,(f)  are 
constants  and  w0  =  u»,.  In  this  case,  a  straight¬ 
forward  but  tedious  calculation  using  equations 
(2),  (11),  and  (24)  through  (28)  gives  y(t)  —  0, 
indicating  complete  cancellation  of  the  primary 
input. 
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2.2  Adaptive  Null  Steering 


We  can  interpret  the  operation  of  the  adaptive 
antenna  aystem  in  terms  of  adaptive  null  steerin'. 
Since  we  are  not  seeking  to  obtain  new  results,  ut 
only  to  interpret  old  ones,  we  simplify  the  mathe¬ 
matics  by  assuming  that  a  >>  1.  Thus,  the  second 
term  of  equation  (32)  is  negligible  compared  to  the 
first,  and  the  third  term  is  negligible  compared  to 
the  fourth.  Decomposing  the  first  and  fourth  terms 
yields 


y(t)  =*  [*(/)  +  J.(0]  ~ 


a  »  1  ,  (35) 

where 

Ji(<)  =  A jit)  cos  (to,f  +  —  0,  +  02\  (36) 

is  a  phase-shifted  version  of  J(f).  The  first  brack¬ 
eted  term  on  the  right  side  of  equation  (35)  can  be 
interpreted  as  the  response  due  to  the  primary 
antenna  with  gains  and  Gpj  in  the  directions  of 
the  desired  signal  and  interference,  respectively. 
The  second  bracketed  term  can  be  interpreted  as 
the  response  due  to  an  equivalent  pattern,  the 
output  of  which  is  subtracted  from  the  primary 
output  to  give  y(f). 

With  this  interpretation,  the  equivalent  can¬ 
cellation  pattern  or  adaptive  beam  has  gains 

G»  ~  |  j  4.  ^  *7».i  <  0  >>  1  •  (37) 

<V  -  (l  +  •  *  »  1  .(38) 

in  the  directions  of  the  desired  signal  and  inter¬ 
ference,  respectively.  The  equivalent  overall  pat¬ 
tern  of  the  adaptive  antenna  system  provides  the 
gains 


w  *  (yGp*  v/c; )  -  ( 1  +  )  Gpt ' 

a  »  1  ,  (39) 


a  »  1  ,  (40) 


in  the  directions  of  the  desired  signal  and  inter¬ 
ference,  respectively.  If  pr  is  small,  equation  (39) 
gives  Gg  =*  G^.  If  pris  so  small  that  prQ 1  «  1, 
equation  (40)  indicates  that  Gy'  is  much  smaller 
than  Gpj.  Thus,  an  approximate  null  can  be  created 
in  the  direction  of  the  interference.  For  this  reason, 
the  action  of  an  adaptive  antenna  system  is  some¬ 
times  called  “null  steering.”  A  pictorial  representa¬ 
tion  is  given  in  figure  3. 


ntm  sjgnm 


Figurir  3.  Adaptive  beam  forming  and  null 
steering. 
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An  analysis  based  upon  the  Wiener-Hopf  equa¬ 
tion,  which  is  derived  in  section  3.3,  shows  that  if 
the  narrowband  condition  given  by  equation  (6b) 
is  not  satisfied,  then  the  performance  of  the  side- 
lobe  cancellei  ?f  figure  1  deteriorate.2  Conse¬ 
quently,  adaptive  systems  intended  to  handle 
wide-bandwidth  signals  must  use  more  elaborate 
processing,  such  as  that  provided  by  delay  lines 
with  multiple  taps,  which  are  discussed  in 
section  3. 

Although  the  major  features  of  the  sidelobe  can¬ 
celler  have  been  determined,  the  theoretical  moti¬ 
vation  for  the  configuration  of  figure  1  has  not  been 
discussed.  For  this  purpose,  and  to  explore  alterna¬ 
tive  configurations,  we  consider  the  classical 
theory  of  adaptive  antenna  systems. 

3.  CLASSICAL  THEORY 

An  adaptive  antenna  system  consists  of  an 
antenna  array,  fixed  processing  elements,  and 
adaptive  elements.  Its  purpose  is  to  remove  ex¬ 
ternally  generated  interference  from  a  receiver 
signal.  The  adaptive  elements  have  one  of  the 
forms  shown  in  figure  4.  In  figure  4(a),  the  weight 
adjustment  mechanism  responds  to  the  output;  in 
figure  4(b),  it  responds  to  the  error,  which  is  the 
difference  between  the  desired  response  and  the 
output.  Ideally,  the  desired  response,  d,  is  the 
received  signal  minus  the  interference  and  noise. 
In  practice,  d  is  a  signal  with  the  general  character¬ 
istics,  but  not  the  detailed  structure,  of  the  signals 
the  antenna  system  is  attempting  to  receive.  The 
forms  in  figure  4  are  not  necessarily  optimal,  but 
are  used  because  of  their  simplicity  and  compat¬ 
ibility  with  computer-controlled  systems.  The 
inputs,  X„  X2,  .  .  .  ,  XN,  are  derived  from  fixed 
processing  elements  and  may  be  either  continuous 
or  discrete-time  signals.  The  inputs  are  applied  to 
weights  W„  W2,  .  .  . ,  W|\|,  which  are  continually 
adjusted  as  a  function  of  the  output  or  the  error. 
The  basic  adaptive  elements  can  be  used  in  a 
variety  of  different  systems. 

i 
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The  inputs  may  be  derived  from  a  tapped  delay 
line,  as  shown  in  figure  5.  Each  of  the  K  antenna 
outputs  is  filtered,  delayed,  and  then  applied  to 
a  line  consisting  of  L  tap  points  and  L  —  1  ideal 
time  delays  of  S  seconds  each.  With  N  =  KL,  the 
adaptive  part  of  figure  5  is  a  special  case  c '  "gure  4.  , 

Each  input  of  figure  4  can  be  the  discrete 
Fourier  transform  at  a  specific  frequency,  say _/j,  of 
an  antenna  output.  With  this  interpretation  for 
figure  4(b),  the  desired  response  is  the  discrete 
Fourier  transform  at  fi  of  the  desired  waveform.  A 
complete  frequency-domain  array-processing  fil¬ 
ter  feeds  K  of  the  adaptive  elements  —  one  for  each 
of  K  discrete  frequencies,  as  illustrated  in  figure  6. 
This  filter  might  be  attractive  if  the  desired  wave¬ 
form  is  a  complicated  function  of  time  but  has  a 
simple  Fourier  transform. 

To  develop  adaptive  algorithms,  we  first  deter¬ 
mine  the  optimal  fixed  weights.  Then,  adaptive 
algorithms  can  be  determined  that  yield  weights 
converging  to  the  optimal  fixed  weights. 

The  derivation  of  the  optimal  fixed  weights 
depends  upon  the  specification  of  a  performance 
criterion  or  estimation  procedure  and  the  modeling 
of  the  signals  received  by  the  antennas.  For  fixed 
weights,  the  optimization  problem  is  illustrated  in 
figure  7.  The  input  and  weight  vectors  are 


X, 

W, 

X, 

W2 

X  = 

• 

,  W  = 

• 

Xn_ 

_Wn_ 

The  input  is  a  function  of  the  desired  signal,  inter¬ 
ference,  and  thermal  noise.  The  scalar  output  is 

y  =  W^X  .  (42) 
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F,gUre  4>  h**ir  *d*pt»e  element*, 


DESIRED 

response 


X 


UNEAR 

y  =  WTX 

FILTER  W 

ESTIMATE  OF 

DESIRED  RESPONSE 

Figure  7.  Optimisation  problea 


where  the  superscript  T  denotes  the  transpose 
operation.  The  output  is  an  estimate  of  the  desired 
response,  d.  The  fiber  or  combiner  of  figure  7  is 
linear  because  the  weights  are  fixed.  Since  their 
characteristics  change  in  response  to  changes  in  the 
characteristics  of  their  inputs,  adaptive  systems 
are  intrinsically  nonlinear.  However,  adaptive  sys¬ 
tems  can  be  designed  that  approximate  the  optimal 
linear  filter  after  a  sufficient  amount  of  adaption. 

Several  different  estimators  of  d  can  be  imple¬ 
mented  by  the  linear  filter  of  figure  7.  If  the 
received  vector,  X,  is  assumed  to  have  a  Gaussian 
distribution,  maximum  likelihood  estimation  of 
the  desired  signal  components  of  the  received 
waveform  leads  to  an  estimator  that  is  a  linear 
function  of  X.  However,  the  Gaussian  assumption 
is  unwarranted  when  interference  is  present.  Esti¬ 
mators  that  depend  only  upon  the  second-order 
moments  of  X  can  be  derived  using  performance 
criteria  based  upon  the  signal-to-noise  ratio  or  the 
mean- square  error  of  the  filter  output.  Conse¬ 
quently,  maximum  likelihood  estimation  is  not 
pursued  further  in  this  paper.  Griffiths 1  compares 
the  maximum  likelihood  filter  and  the  Wiener 
filter,  which  is  derived  below. 


3.1  Signal-to-Noiac  Ratio  Criterion 

A  reasonable  design  criterion  for  adaptive  radar 
systems  is  the  maximization  of  the  probability  of 
detection  for  fixed  probability  of  false  alarm.4  This 
maximization  has  been  shown  to  be  equivalent  to  a 
maximization  of  a  generalized  signal-to-noise  ratio, 


which  ia  itself  an  intuitively  appealing  design 
criterion. 

Consider  the  maximisation  of  the  signal-to- 
noise  ratio  at  the  output  of  the  adaptive  element  of 
figure  4a.s  By  noise,  we  mean  the  combined 
thermal  noise  and  externally  generated  interfer¬ 
ence.  We  assume  that  the  input  Vector,  X,  ia  the 
sum  of  two  components. 


X  -  H  d  +  n  , 
where  H  is  a  vector  of  constants. 


(43) 


n , 

nt 


(44) 


"N 


is  a  vector  of  noise  components  that  are  derived 
from  stationary  processes,  and  d  is  the  determin¬ 
istic  desired  signal.  We  treat  all  vector  components 
and  scftiars  as  complex  variables.  Consequently, 
our  results  apply  to  adaptive  elements  fed  by  the 
system  of  figure  6  and  to  narrowband  systems  for 
which  d  is  the  complex  envelope  representation  of 
the  desired  signal. 

The  component  of  the  output,  y,  that  is  due  to 
the  desired  signal  is 


L  l  Griffiths,  A  of  fduhiiimentimol  Wiener  and  Maximum- 

likelihood  Filter*  for  Amtenno  Arroyo,  Proc.  IEEE  (letterob  55,  (November 
19671*04$. 

4 
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yt  -  WTH  A  , 


(45) 


&  P.  Appttbmum,  Adoptive  Arroyo,  IEEE  Tram*.  Anunnao  Prop.,  AP-2* 
(September  1976b  «S 
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and  the  component  of  y  that  is  due  to  the  noise  is 
y«  “  WTn  »  nTW  .  (46) 

At  a  single  instant  of  time,  the  signal  power  at 
the  output  is 

R*  -  W1  »  IwTHdl*  .  (47) 

Let  a  prime  denote  the  conjugate  transpose  opera¬ 
tion  and  an  asterisk  denote  the  conjugate  opera¬ 
tion.  The  mean  noise  power  in  the  output  is 

R„  -  E[(nTW)'(nTw)]  -  W'RnnW  , 

(48) 

where  E(x]  is  the  expected  value  of  x  and 
Rnn  *■  E[n*nT]  is  the  correlation  matrix  of  the 
noise.  The  signal-to-noise  ratio  is 

p  =  (WTH),(WTH)|dl* 

Rn  W'RnnW 

(49) 

Since  R„  >  0  in  practical  cases,  Rnn  is  positive 
definite.  From  its  definition,  Rnn  ia  Hermitian. 

A  positive  definite  Hermitian  matrix  has  posi¬ 
tive  eigenvalues.  Thus,  it  can  be  diagonalixed  to  the 
identity  matrix.  We  define  the  nonsingular  N  X  N 
matrix, 

e,  e,  ®n 

D  =  —  ...  ,  (50) 

v^A,  V  A,  v/An 


where  the  A,-  are  the  eigenvalues  and  the  c,  are 
the  corresponding  orthonormal  eigenvectors. 
Straightforward  matrix  algebra  yields 

I  =  D'  Rnn  D  ,  (51) 

where  I  is  the  identity  matrix. 


Combining  this  result  with  equation  (48)  gives 
Rn  -  W'  (D')-‘  D'  Rnn  DD‘  W  - 

(D-‘w)'(D-,W)  -  II  D-‘W  II 1  ,  (52) 

where  ||  x  ||  ■*  x’  x  is  the  Euclidean  norm  of  the 
veca>i  x.  Using  the  Cauchy-Schwarz  inequality, 
equation  (47)  yields 

R,  -  WT(DT)  'DTH  d[  - 

j(D-,W)1OTw|*  - 

|  (D"'W)'D'H*  fUl*  ^ 

IId-'wI*  II  D'H*  II1  d|*  .  (53) 

Equality  in  this  equation  is  attained  if 

D  'W  -  cD'H*  ,  (54) 

where  c  is  an  arbitrary  constant.  Substituting 
equations  (52)  and  (53)  into  equation  (491,  it 
follows  that 

P  <  IID'H*  IP  I  d\*  ,  (55) 

The  maximum  value  of  P  is  attained  if  equation 
(54)  is  satisfied.  Thus,  the  optimal  choice  of  W  is 

W,  -  cDD'H*  .  (56) 

Using  equation  (51),  this  relation  becomes 

W,  -  cRnnH*  .  (57) 

The  value  of  P  corresponding  to  optimal  weights  is 

p.  -  iiDu-irur  - 

H7toD'H*|i/r  -  mr^ldr  .  (58) 

In  order  to  implement  equation  (57),  we  must 
know  the  values  of  the  components  of  H.  Since 


adaptive  systems  have  more  then  one  antenna  (or 
array  element),  this  knowledge  implies  that  the 
direction  of  arrival  of  the  desired  signal  is  known. 

3.2  Mean-Square  Error  Criterion 

The  mobt  widely  used  method  of  estimation  is 
based  on  the  minimization  of  the  mean-square 
error.*  Instead  of  assuming  a  known  direction  of 
arrival,  we  assume  that  a  facsimile  of  the  desired 
vesponse  is  supplied.  We  assune  that  X  and  d  are 
derived  from  stationary  stochastic  processes. 
Similar  results  can  be  obtained  by  modeling  X  and 
d  as  continuous-time  or  discrete-time  determin¬ 
istic  variables  and  minimizing  the  time-  or  sample- 
averaged  square  error. 

Because  it  is  convenient  and  because  the  results 
are  usually  applied  to  systems  of  the  form  of  figure 
5,  we  assume  that  all  scalars  and  matrix  com¬ 
ponents  are  real.  Referring  to  figure  4(b),  the 
difference  between  the  desired  response  and  the 
output  is  the  error  signal, 

e  =  d  -  WTX  .  (59) 

We  sqt  are  both  sides  of  this  equation  and  take  the 
expected  value.  Regarding  W  as  fixed,  and  noting 
that  WTX  =  XTW,  the  mean-square  error  is 

Eje2]  =  E[(rf  -  WrX)2]  = 

E  M  +  WTR„W  -  2WTRw  ,  (60) 

where 

Rxx  -  e[xxt]  (61) 

is  the  symmetric  correlation  matrix  of  X,  and 

R*d  =  E[Xd]  (62) 

is  a  vector  of  cross  correlations  between  the  input 
signals  and  the  desired  response.  Taking  the 
- 
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gradient  of  equation  (60)  with  respect  to  W  yields 

7  Etc1]  =  2BXXW  -  2Rxd  .  (63) 

Setting  the  gradient  equal  to  zero,  we  obtain  a 
necessary  condition  for  the  optimal  weight  vector 
W„.  In  practice,  R**  is  positive  definite.  Thus,  it  is 
nonsingular  so  that  we  obtain 

W0  =R- R^  .  (64) 

This  equation  is  the  Wiener-Hopf  equation  for  the 
optimal  v/eight  vector.  The  associated  linear  filter 
is  called  the  Wiener  filter.  To  show  that  W0 
produces  the  minimum  mean-square  error,  we  set 
W  =  W0  in  equation  (60)  to  obtain  the  mean 
square  error  corresponding  tc  W0* 

Er  =  E[</*J  —  WjR^Wo  - 

E[d*]  ~  •  (65) 

Substituting  equation  (65)  into  equation  (60)  and 
using  equation  (64)  to  eliminate  R*,/,  we  may 
express  the  mean-square  error  in  the  form. 

EM=  Er  4-  (W  —  W0)T  R**(W  —  W0)  .  (66) 

Sipice  R**  is  positive  definite,  this  expressio, 
shows  that  W0  is  a  unique  optimal  weight  vector 
and  Er  is  the  minimum  mean-square  error.  Equa¬ 
tion  (66)  describes  a  multi-dimensional  quadratic 
function  of  the  weights  that  can  be  visualized  in  two 
dimensions  as  a  bowl-shaped  surface.  The  purpose 
of  adaptation  is  tc  continually  seek  the  bottom  of 
the  bowl. 

3.3  Steepest  Descent 

In  the  implementation  of  either  equation  (57)  or 
equation  (64),  the  presence  of  interference  means 
that  the  correlation  matrices,  R-^n  or  R^,  are 
unknown.  For  this  reason,  approximations  cf  the 
equations  must  be  determined.  Since  the  computa¬ 
tional  difficulties  of  the  matrix  inversions  are 
considerable  when  the  number  of  weights  is  large 
or  time-varying  signal  statistics  require  frequent 
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computations,  most  of  the  approximations  devel¬ 
oped  involve  recursion  relations  and  gradient  or 
random  search  techniques.7  However,  for  sampled 
data  systems,  direct  implementation  of  the  equa- 
tii  as  after  suitable  approximations  of  the  matrix 
elements  has  been  shown  to  exhibit  rapid  con¬ 
vergence  in  certain  circumstances.* 

Because  the  theory  has  been  extensively  inves¬ 
tigated  and  it  has  often  been  implemented,  we 
consider  the  approximation  of  equation  (64)  by  the 
method  of  steepest  descent.  The  Approximation  of 
equation  (37)  implies  a  similar  adaptive  weight 
adjustment  mechanism.*  We  shall  treat  X,  W,  and 
d  as  discrete-time,  sampled-value  variables  and 
use  the  index  k  to  denote  a  particular  sampling 
instant  or  adaptation  cycle.  The  results  for  contin¬ 
uous-time  systems  are  analogous. 

In  the  method  of  steepest  descent,  the  weight 
vector  is  changed  along  the  direction  of  the  nega¬ 
tive  gradient  of  the  mean-square  error.  Using 
equation  (63),  we  obtain 

W(/fc  +  l)  =  W(*)-MVE[eJ]  = 

W(*)-2M[Rr;tW(*)-  R^]  ,  (67) 

where  the  scalar  constant  /x  controls  the  rate  of 
convergence  and  stability.  The  adaptation  cycle 
begins  with  an  arbitrary  initial  weight.  As  it  stands, 
this  method  eliminates  the  need  to  compute  the 
inverse  of  R**,  but  still  requires  approximation  of 
R**  and  R^. 

The  Griffiths  algorithm9  result*  if  R**  is  re¬ 
placed  by  X(A)XT(k)  and  is  directly  estimated, 
that  is,  we  use 

W((k+1)  =  W(A)-2M[X(*)XT(A)W(*)-R*rf] 


L  E.  Brennan  and  i  S.  Heed,  t'hcor>  of  Adaptive  Hadar ,  /fH  Trans. 
Aerospace  Electron.  Syst.,  AES-9  (March  I 973),  237. 

ft  W id  row  and  J.  M.  McCool,  A  Comparison  of  Adoptive  Algorithms  Based  on 
Oh  Methods  of  Steepest  Descent  and  fandom  Scorch,  IE  EE  Trans.  Antennas  Prop., 
A  P-24  (September  1976),  6IS. 

1.  S.  Heed  et  aL  Rapid  Convergence  Rates  in  Adaptive  Arrays,  IEEE  Trans. 
Aerospace  Electron.  Syst .  I ES-IO  (November  1974),  853. 
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L  J.  Griffiths,  •!  Simple  Adaptive  Algorithm  for  Real-Time  Processing  in 
.4  *  fen  *  a  -Imivi,  Pmc.  IEEE,  57  (October  196%  Ib9b 


-  W(A)  —  2fi[\(k)y(k)  —  fk*,j  ,  (68) 

where  is  an  estimate  of  R,</.  The  problem 
with  this  algorithm  is  that  an  accurate  estimate  of 
Rxrfis  difficult  if  the  amplitude  of  the  desired 
signal  is  unknown. 

The  Widrow-Hoff  algorithm,4  also  known  as  the 
least-mean-square  (LMS)  algorithm,  approxi¬ 
mates  R^  by  X(A)d(A)  so  that 

W (A  +  1)  =  W (k)  -  2p[X(k)y(k)  -  X(A)d(A)] 

=  W(A)  -  2fi£(k)X(k)  .  (69) 

"'his  iteration  rule  states  that  the  next  weight 
vector  is  obtained  by  adding  to  the  present  weight 
vector  the  input  vector  scaled  by  the  value  of  the 
error.  It  can  be  shown  that  for  a  sufficiently  small 
value  of  ju,  the  mean  of  the  weight  vector  converges 
to  the  optimal  value  given  by  the  Wiener-Hopf 
equation.  Figure  8  shows  a  block-diagram  repre¬ 
sentation  of  equation  (69)  for  one  component  of  the 
weight  vector. 

In  analog  implementations  of  continuous  sys¬ 
tems,  equation  (69)  becomes  the  differential 
equation, 

-^W (I)  =  -  2 ps(t)X(t)  .  (70) 

Equivalently,  if  W(0)  =  0  , 

W(f)  =  —2p  £  e(r)X(r)dr  .  (71) 

Figure  9  shows  a  realization  of  one  component  of 
this  er  >ation. 

One  potential  problem  with  the  Widrow-Hoff 
algorithm  can  be  explained  in  terms  of  the  overall 
equivalent  radiation  pattern  of  the  antenna  array. 
Adaptation  causes  the  array  to  form  a  beam  in  the 
direction  of  the  desired  signal,  while  reducing  or 
nulling  the  array  response  in  the  direction  of  the 
interference.  If  the  desired  signal  is  not  transmitted, 
the  adaptation  may  proceed  to  create  a  reduced 
response  in  the  direction  of  this  signal  so  that 

*8.  WiJrownmL  AlmptivrAiUmno System*.  Proe.  IEEE,  55  {Oeeember  IViT), 
21*3. 
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Figure  8.  Digital  realisation  of  Widrow-Hoff  algorithm. 


Figure  9.  Analog  realisation  of  Widrow-Hoff  algorithm. 


communication  is  hindered  or  disabled  when 
transmission  is  resumed.  To  alleviate  this  prob¬ 
lem,  a  pilot  signal  may  be  used. 

A  pilot  signal  is  a  vector,  each  component  of 
which  is  an  estimate  of  the  desired  signal  that 
would  appear  at  an  antenna  output  if  no  interfer¬ 
ence  or  noise  were  present.  In  other  words,  the  pilot 
signal  is  constructed  to  have  spectral  and  direction¬ 
al  characteristics  similar  to  those  of  the  incoming 
signal  of  interest.  Two  different  algorithms  that  use 
the  pilot  signal  have  been  proposed.  In  the  two¬ 
mode  algorithm,  adaptation  alternates  between 
adapting  only  on  the  pilot  signal  and  adapting  only 
on  the  signals  actually  received  by  the  array.  The 
one-mode  algorithm  adapts  on  the  sum  of  the  pilot 
signal  and  the  received  signal.  Details  of  the 
performance  of  the  pilot  signal  method  can  be 
found  in  the  literature.*-  * 

*&  WOn mHal  Adoptim  Anmno  Srtirmt,  Proe  /££  SS  fOretmb rr  1967), 
}l*r 

*L  J.  (iriffitki.  ,4  Simple  Adoptior  Algorithm  for  ftrobTimr  Procruing  in 
Antmn*  Arrmyt,  Proc.  IEEE,  S7  (Ortobtr  1969),  1696, 


The  desired  response,  d,  is  required  at  each 
iteration  of  equation  (69).  An  estimate  of  the 
desired  response,  denoted  by  <J,  can  be  generated 
for  some  types  of  communications  by  feedback 
systems  of  the  forms  shown  in  figure  10.  For 
effective  operation,  et  does  not  have  to  be  a  perfect 
replica  of  the  desired  response;  if  it  were,  we  would 
not  need  the  adaptive  antenna  system.  However, 
the  Wiener-Hopf  equation,  which  we  would  like  to 
approximately  satisfy,  indicates  that  i  should  be 
such  that 


R*a  ^  R rrf  (72) 

Stable  operation  is  achieved  only  if  &  is  inde¬ 
pendent  of  the  amplitude  of  the  component  propor¬ 
tional  to  the  desired  signal  in  the  array  output.  If  it 
is  not,  the  array  weights  may  increase  indefinitely 
or  decrease  to  zero.  The  presence  of  the  limiter  in 
figure  10(b)  controls  the  amplitude  of  H.  The  delay 
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Figure  10.  Generator*  of  desired  mponie. 

duj  lo  the  feedback  generation  of  <2  should  be 
considerably  less  than  the  data-modulation  period. 
Details  of  the  operation  of  an  adaptive  array  in  a 
spread-spectrum  communication  system  and  with 
feedback  generation  of  the  desired  response  are 
described  by  Compton. For  this  system,  the 
interference  suppress.on  due  to  the  adaptive  array 
is  available  in  addition  to  the  suppression  afforded 
by  the  waveform  processing. 

4.  ADAPTIVE  NOISE  CANCELLING 

The  adaptive  element  of  a  noise  cancelling 
system1 1  is  a  opecial  case  of  figure  i{»).  However,  a 
noise  canceling  system  has  such  distinctive 


ft  J.  Compton,  An  Ainptipe  Army  tn  •  SprmdSpectrum  Communiction 
System,  Prvc.  i  'EE,  66  (Msrrh  1970k  289. 

B.  Widrvw  et  ni  Adoptive  Noise  Cnneellinp:  Principles  snd  Applications, 
Proc.  IEEE,  63  ( December  J975J.  1692. 


features  that  it  is  treated  separately  in  this  section. 
In  adaptive  noise -cancelling  systems,  illustrated  in 
figure  1 1,  the  desired  response  is  replaced  by  the 
primary  input,  which  is  derived  from  a  separate 
antenna.  The  primary  input  is  excluded  from 
weight  vector  calculations.  Thus,  the  optimal 
weight  vector  for  minimizing  the  mean  square 
error  is  given  by  the  Wiener-Hopf  equation, 

V0  —  Ri*  .  (73) 

where  we  interpret  as  the  correlation  matrix  of 
the  reference  inputs,  and  R*p  is  the  correlation 
vector  of  the  reference  inputs  with  the  primary 
input. 

The  system  output  for  the  adaptive  noise¬ 
cancelling  system  is  identical  to  the  error  signal. 
We  show  that  minimizing  the  mean-square  output 
is  approximately  equivalent  to  causing  the  output 


to  be  u  ir  initnum  mean-square  error  estimate  of  the 
desired  signal,  d,  if  the  adaptive  filter  output,  y, 
is  such  that  it  has  negligible  correlation  *ith  d.  We 
assume  that  the  primary  input  is  the  sum  of  the 
desired  signal  and  uncorrelated  noise,  n*  Refer¬ 
ring  to  figure  1 1,  the  canceller  output  is 

c  -  d  +  n,  -  y  .  (74) 


We  assume  that  all  signals  are  derived  from 
stationary  stochastic  processes.  Taking  the  ex¬ 
pected  value  of  the  square  of  both  sides  of  equation 
(74)  gives 

Efc]  =  E[o”]  +  E[(n.-y)lJ  + 

2E[dn,]  -  2E[dy]  .  (75) 
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If  J  ia  uncorrelated  with  n«  and  haa  negligible  cor¬ 
relation  with  y,  thin  equation  reducee  to 

£[•>]- E[*]  +E[(»,-,)‘]  (76) 

The  aignal  power,  E[rf*],  ia  unaffected  aa  the  adap¬ 
tive  filter  ia  adjusted  to  minimite  E[e*].  Conse¬ 
quently,  E[(n,  -  y)1]  ia  minimised  when  E[eJ] 
ia  minimised.  From  equation  (74),  it  follows  that 
E[(e  —  </)’]  is  also  minimised.  We  conclude  that 
the  adaptive  noise  canceller  output  is  a  minimum 
mean-square  error  estimate  of  the  desired  signal. 

Since  the  Wiener-  Hopf  equation  is  valid,  the 
Widrow-Hoff  algorithm  applies.  Consider  a 
continuous-time  system  with  a  two-dimensional 
reference  vector,  XT  =  [X„  X,],  where  X2(t)  is 
the  result  of  passing  X|(f)  through  a  quarter- 
wavelength  delay.  Equation  (71)  then  implies 
the  implementation  of  Figure  1.  Thus,  we  have 
shown  how  this  configuration  arises  and  why  its 
performance  can  be  expected  to  be  nearly  optimal. 
Using  the  Wiener-Hopf  equation,  Widrow"  derived 
a  Z-transform,  sampled-data  version  of  equation 
(34). 


such  as  removing  periodic  interference  from 
spread-spectrum  communications.  The  adaptive 
notch  filter  has  the  form  depicted  in  figure  12, 
which  results  from  figure  1 1  when  the  reference 
inputs  an'  derived  from  a  tapped  delay  line  that  is 
fed  by  the  delayed  primary  input.  The  delay  A  is 
sufficiently  long  to  cause  the  desired  wideband 
signal  component  of  X,  to  become  uncorTelated 
with  the  corresponding  component  of  the  primary 
input.  However,  because  of  their  nature,  periodic 
interference  components  in  X,  usually  retain  cor¬ 
relation  with  the  corresponding  components  of  the 
primary  input.  As  a  result,  the  tapped  delay  line 
adaptively  forms  a  filter  such  that  y  is  composed 
primarily  of  the  periodic  interference  components. 
If  the  filter  bandwidth  is  sufficiently  narrow,  the 
system  output  is  the  desired  wideband  signal  with 
little  distortion. 


As  shown  in  section  1 .  adaptive  noise  cancelling 
results  in  a  small  amount  of  cancellation  of  the 
desired  signal.  In  terms  of  the  overall  antenna 
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Figure  1 2.  Adaptive  notch  filler. 


The  adaptive  noise  canceller  can  be  used  as  an 
adaptive  notch  filter,  which  is  useful  in  applications 


Widrow  tt  *4  Adaptive  .Voile  CanrrtUnp:  Principirt  and  Applications, 
* or  IEEE,  M.  (DK*m6*r  197$k  926. 


pattern,  this  cancellation  is  caused  by  the  changes 
in  sensitivity  of  the  main  beam  in  the  steering 
direction.  These  changes  can  occur  because  the 
mainlobe  pattern  is  not  constrained  by  the  adap¬ 
tive  process. 


23 


5.  CONSTRAINED  MINIMUM 
POWER  CRITERION 

A  performance  criterion  that  inherently  limits 
the  cancellation  of  the  desired  signal  is  the  con¬ 
strained  minimum  power  criterion,  which  is  ap¬ 
plied  to  systems  with  adaptive  elements  of  the  form 
of  figure  4(a).  We  assume  that  the  components 
of  X  are  derived  from  stationary  stochastic 
processes.  The  cross-correlation  matrix  is  defined 
by  R*,  ™  E[X*XT].  The  criterion  requires  the 
minimization  of  the  mean  output  power, 


For  the  system  of  figure  5,  the  constraint  is 
somewhat  different.  We  assume  that  identical 
desired  signal  components  arriving  on  a  plane 
wavefront  appear  at  the  first  K  taps  simultan¬ 
eously.  Consequently,  the  filter  response  to  the 
signal  is  equivalent  to  the  response  of  a  single 
tapped  delay  line  in  which  each  weight  is  equal  to 
the  sum  of  the  weights  in  the  corresponding 
vertical  column  of  the  adaptive  element.  These 
summation  weights  in  the  equivalent  line  are 
seated  to  give  a  desired  frequency  response. 
Thus,  we  -squire  that 


Elly:2]  -  W'l^W  ,  (77) 

subject  to  the  constraint 

C'W«F  ,  (78) 

where  F  is  a  specified  constraint  vector  and  C  is  a 
matrix  associated  with  the  constraint.  The  specifi¬ 
cation  of  C  and  F  depends  upon  the  form  of  the 
adaptive  system. 


C;  W  -  /  .  :  -  1,  2,  ..  .  .  ,  L  ,  (80) 

where  /•  is  a  desired  summation  weight  and  C|  is 
an  N-dimensional  vector  of  the  form 

Cf  -  [0  0  ...  1  1  ...  1  0  0  ...  0]  .(81) 

The  l's  correspond  to  the  ah  vertical  column  of  K 
weights.  To  put  equation  (80)  in  matrix  form,  the 
constraint  matrix  is  defined  as 


For  the  system  of  figure  6,  equations  (77)  and 
(78)  are  applied  to  each  summer  output.  Consider 
a  specific  summer  and  let  W  represent  the  associ¬ 
ated  complex  weight  vector  and  X  the  associated 
complex  input  vector.  Assuming  that  the  steering 
delays  are  such  that  the  antenna  array  pattern 
points  in  the  direction  of  the  desired  signal,  the 
elements  of  X  differ  from  each  other  only  because 
of  the  interference  and  noise.  Thus,  if  we  require 
equation  (78)  with 

C'=  [l  1  ...  l].F  =/,  (79) 

where/ is  a  scalar,  this  constraint  forces  the  de¬ 
sired  signal  component  of  the  summer  output  to  be 
a  conventional  coherent  sum.  Minimization  of  the 
output  power  then  minimizes  the  interference 
power  obtained  from  directions  other  than  the 
main  beam  direction.  Note  that  W  has  2N  degrees 
of  freedom  because  of  the  N  magnitudes  and 
N  phases  of  its  elements.  The  complex  constraint 
equation  removes  two  of  these  degrees  of  freedom, 
but  leaves  2N  —  2  to  be  used  in  limiting  the 
interference. 


C  =  [C,  C,  ...  CL]  ,  (82) 

which  has  the  dimensions  N  by  L  TheL-dimen- 
sional  vector  of  weights  of  the  equivalent  line  is 


(83) 


A 


With  these  definitions,  the  constraints  can  be 
combined  into  a  single  matrix  equation  of  the  form 
of  equation  (78).  Once  again,  the  constrained 
minimization  ensures  that  the  output  power  from 
all  directions  except  boresight  is  minimized. 

Having  shown  how  the  constraint  equation  is 
constructed,  we  proceed  to  a  derivation  of  the 
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where 


optimal  weight  vector,  aaauming  that  the  weighta 
are  fixed  constant*.  For  simpl;  Sty,  we  aaaume  that 
all  acalar*  and  matrix  component*  are  real.  By  the 
method  of  Lagrange  multipliers,  we  minimise 

H  -  WTR*,W  +  XT(CTW  -  F)  ,(84) 

where  X  ia  the  vector  of  Lagrange  multipliers. 

Taking  the  gradient  of  equation  (84)  with  re¬ 
spect  to  the  vector  W  yields 

VH  -  2R*,W  +  C  X  .  (85) 

A  necessary  condition  for  the  minimum  is  deter¬ 
mined  by  setting  VH  equal  to  sero.  Thus, 

w.  —  ~  ~  R*i  CX  (86) 

Since  W,  must  satisfy  the  constraint,  we  substi¬ 
tute  equation  (86)  into  equation  (78)  with  the 
result, 

-  j  CX-  F  (87) 

This  equation  can  be  solved  for  X.  Substituting 
the  solution  into  equation  (86)  gives 

W,  -  Ht,  C(CTR^C)-'  F  ,  (88) 

where  we  assume  that  the  indicated  inverse  exists. 
Frost’s  adaptive  algorithm,12  which  is  based  upon 
approximating  equation  (88),  is 

W<0)  -  B  , 

W(*  +  1)  - 

A[W(*)  -  ny(k)\(k,]  +  B  ,  (89) 


O.L.  Fruit,  -In  AUonthm  for  hnroriy  Co*strmin+d  idoptii*  -frrtn  /Vocrj- 
un*,  /Vor  IEEE,  bO  1 4ugu*t  1972 A  926 


A  -  I  -  C(CTC)-‘  CT  , 

B  -  C(CTC)-'  F  .  (90) 

Hie  algorithm  ensures  the  satisfaction  of  the  con¬ 
straint  after  each  iteration  and  has  other  desirable 
properties.  Neither  pilot  signals  nor  generators  of 
the  desired  response  are  needed  in  constrained- 
minimum-power  systems. 

In  the  development  of  the  adaptive  algorithms, 
it  has  always  been  assumed  that  the  noise  proc¬ 
esses  are  stationary.  However,  the  main  purpose  of 
adaptive  systems  is  to  automatically  adjust  to 
nonstationary  inputs,  especially  when  the  sto¬ 
chastic  properties  are  unknown  a  priori.  For¬ 
tunately,  experimental  results  and  tentative  theo¬ 
retical  results  seem  to  indicate  that  adaptive  algo¬ 
rithms  retain  their  usefulness  in  many  realistic 
environments  with  nonstationary  inputs. 11 

In  addition  to  the  stochastic  properties  of  the 
inputs,  we  have  tacitly  assumed  plane-wave  sig¬ 
nals,  an  ideal  propagation  medium,  and  distortion¬ 
less  receivers.  Since  these  assumptions  may  not  be 
valid  under  actual  operating  conditions,  Vural  has 
investigated  the  effects  of  signal,  medium,  and 
system  deviations  from  the  usual  assumptions.14 


Adaptive  antenna  systems  have  many  applica¬ 
tions.  In  general,  adaptation  is  potentially  helpful 
for  communications  or  radar  when  the  desired 
signd  and  the  interferences  are  distinguishable  a 
priori .  The  discriminant  may  be  the  direction  of 
arrival,  a  waveform  characteristic,  or  even  signal 
polarization.1 


i 
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